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Introduction 


The zodiacal cloud is populated by cometary, asteroidal (Leinert and Grim, 1985) and recently, 
interstellar, dust particles (Griin et al., 1994) each with its own distinct particle properties and 
spatial distribution. One of the main limitations in understanding the structure and the origin of 
the cloud is that all observations yield the integrated flux along any given line of sight, making 
it difficult to identify both the type of source and the location of the radiation. This problem 
can be overcome by studying individual structures in the cloud that are associated with a 
specific source, for example, using the zodiacal dust bands that have a asteroidal origin can 
help us to model the asteroidal component of the cloud. 

This reports presents the latest results from the ongoing project to study and understand the 
fine structure of the zodiacal cloud. In this project we have concentrated upon the near-infrared 
(near-IR) observations of the zodiacal cloud by the Diffuse Background Infrared Experiment 
(DIRBE) aboard the Cosmic Background Explorer satellite (COBE). At thermal wavelengths 
(> 4pm) the radiation from the zodiacal cloud is the dominant source of the background flux 
brought about by the heat radiated from the solar system dust particles. However, in the near- 
IR wavebands the signal is generated by scattered sun light and the flux decreases by more than 
an order of magnitude from the thermal emission. Therefore, it is a challenge to study the 
structures, easily detected in thermal emission, in the near-IR. 

DIRBE-COBE Data Set in the near-IR bands: 

The DIRBE instrument aboard the COBE satellite was a ten band photometer measuring 
absolute fluxes by chopping between the sky and an internal reference source. The short term 
and long term photometric stability has been reported to be better than 1% (Spiesman et al., 
1995). The detectors have a 0.7° by 0.7° field of view, and the satellite covered solar 
elongation angles between 64° and 124°. The instrument was in operation for 10 months from 
December 1990 and the initial data products were released in 1994 (COBE Explanatory 
Supplement). When the cryogen ran out, data collection from six of the 10 bands was 
discontinued. Fortunately, the instrument continued to map the sky in the four near infrared 
spectral wavebands at 1.25, 2.2, 3.5 and 4.9 microns. At the end of 1993 when the mission 
drew to a close, the instrument had acquired a vast quantity of data covering three years. These 
data have a lower signal to noise ratio because a) the contribution of stellar point sources 
emission that negligible at longer wavelengths but cannot be overlooked in the spectral bands, 
b) the decrease in the zodiacal brightness by at least an order of magnitude and c) the increase 
in the ambient temperature after the depletion of the cryogen. However, this is a unique data set 
consisting of a redundantly sampled map of the zodiacal sky because DIRBE covers nearly half 
the sky in one day. The same source, Sirius, was used for absolute calibration for all the bands 
which can give a higher accuracy in color determination. The data set used is the weekly 
averaged data that has 41 weeks of redundant coverage of the sky. 

The first step in the analysis of the zodiacal cloud is to simulate scans in ecliptic latitude at 
constant solar elongation angle. As the line-of-sight of the telescope passes closer to the sun 
with decreasing solar elongation angle, the overall flux of the background increases as we 
sample hotter particles along the line of sight. Therefore, a set of scans each at a constant solar 
elongation angle samples particles at different radial distance from the sun at least near the 
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ecliptic plane. Figure. 1 shows a sample scan created at 1.2pm. Although the zodiacal cloud is 
much weaker the underlying shape is evident. In addition there are numerous point sources at 
all latitudes and the galactic plane is the broad peak below the ecliptic. 


Total Signal 



Figure It A sample scan at a constant solar elongation of 75° in band 1 (1.2 microns). 


The second step in the analysis is to separate the smooth low frequency component of the 
zodiacal cloud from the high frequency fine structure. We use a fourier filter technique which 
separates out the low frequency background using a coefficient between 0 and 1 that 
determines the number of frequency components that are included. The filter also incorporates 
a ‘parzen’ windowing method that decreases frequency leakage. Before the filter is applied to 
the scan, the signal is pre-processed to remove the very bright point sources from the image. 
This is important for studying the high frequency component. Figure 2 and 3 show both 
components obtained after applying the median filter and the parzen Fourier filter to the scan in 
Figure 1. 

Initially we had hoped to apply SKY model by Cohen to estimate the sky contribution to the 
the stellar sources but the model is applicable only in band 3 (at 3.4 microns) and is a very 
strong function of the magnitude limit chosen for the stars. 
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Systematic Variations of the Zodiacal Light: 

We first present the results of the study of the smooth background flux in scattered light. The 
data analysis was done in the following steps. 

a) Each week of data was first divided into two sets identified as Leading and Trailing data. 
As the COBE satellite was in a polar orbit, half the observations were taken in the 
descending orbit (Leading) while the other half in the ascending (Trailing) orbit. It is 
important to distinguish the two observations due to the presence of the asymmetry caused 
by the Earth’s ring. 

b) The scans were simulated at elongation angles between 65° and 124° at every 1° for each 
week of data. 

c) The ‘parzen’ filter (described in the previous section) was applied to all the scans and peak 
of the smooth background was calculated for each week of data. The weeks for which the 
Galactic plane fell very near the ecliptic were not used. 

d) For each week, the peak flux (near the ecliptic) was plotted as a function of the solar 
elongation and the a normalized value of the peak flux at a constant solar elongation of 90° 
was calculated, again using a simple first of second order polynomial fit. A value of 90° is 
chosen as it has some special symmetry characteristics. The value of ecliptic longitude also 
remains along the scan at this elongation. In addition the line-of-sight 1 80° away in ecliptic 
longitude is also at 90° solar elongation. 

e) The annual systematic variation of the zodiacal background is measured by plotting the 
peak flux as a function of time. 

Figure 3, 4 and 5 show the results of the extensive data reduction by the variation of the peak 
zodiacal flux in the three near-IR wavebands as function of ecliptic longitude of the Earth 
which corresponds to time of year. The filled circles denote data taken in the leading direction 
(in the direction of the Earth’s motion along its orbit) while the open triangles represent trailing 
data. We see that there is a systematic variation of the flux which is due to the Earth moving 
through a rotationally asymmetric zodiacal cloud in an elliptical orbit. The flux is a minimum 
close to the perihelion of the Earth at 104°. 

The other important effect that was being sought was the increased flux in the trailing direction 
at all times of the year. This is caused by a cloud of dust particles that follow the Earth in its 
orbit. The dust particles are a part of the Earth’s circumsolar ring formed by resonant trapping 
of the low eccentricity asteroidal dust particles as they drift towards the sun from the Asteroid 
Belt due to P-R drag. We see that at 1 .2 and 2.2 microns there is no evidence of the ring while 
there is a slight indication of the asymmetry which is not quantifiable within the errors. 
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Figure 3,4, Si The annual variation of the zodiacal peak flux in three DIBE near-IR bands. 


Zodiacal Dust Bands: 

The zodiacal dust bands are high frequency structures clearly visible as symmetrically placed 
‘bumps’ above the smooth zodiacal background and it is conventional to study them 
separately. The dust bands are formed by the collisional evolution of asteroids in Hirayama 
families (Dermott et al., 1984). Asteroids within a family have similar semi -major axes, 
eccentricities and inclinations because they are formed by the breakup of a single parent 
asteroid. 

Although the asteroids themselves are located in the main belt, the dust produced from their 
collisional evolution extends all the way to the sun due to solar radiation pressure forces in 
time scales similar to their orbital precession. This leads to the formation of a three 
dimensional torus, the edges of which have a higher density of orbits and rise above the 
background zodiacal cloud to form two bands, above and below the ecliptic, for each family. 
Thus, the visible part of the dust bands is only 10% of the total contribution from the families 
(Dermott et al., 1994) as most of the dust is indistinguishable from the background. The dust 
bands are extremely important features because we can identify them with a specific family, 
which helps us to isolate and study the not only the asteroidal component of the zodiacal cloud 
but also identify characteristics of dust particles belonging to individual families. 
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Figure 6: The dust bands in four infrared wavebands of IRAS with resolution of 2 arcminutes. The 
zodiacal cloud has been filtered out using a fast Fourier filter. 



Figure 7: A scan from the DIRBE data in the 25 micron waveband after filtering the smooth zodiacal 
cloud shows evidence of the two pairs of dust bands present in the near-ecliptic band pair. 


There are two prominent band pairs, one at ±10° and the other near the ecliptic (Figure 6). The 
method of co-adding is so effective that the near-ecliptic bands have been further resolved into 
two bands at 0.01% of the total background signal. The two ecliptic band pairs are associated 
with two distinct families, Koronis and Themis (Sykes, 1992), while the ten-degree band is 
presumed to be associated with Eos although conclusive evidence has yet to be found. 
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(Dennott et al., 1992). Co-adding also becomes invaluable for extracting dust bands from the 
100 micron waveband where the signal from the Galaxy interferes strongly with the solar 
system emission and it is impossible to detect the dust bands in an individual profile or even 
after co-adding over a few weeks of data. A sample of the low resolution COBE dust bands is 
shown in Figure 7 with the location of the primary asteroid families near the Ecliptic. 

Dust Bands in Scattered Light: 

The zodiacal dust bands in scattered light have been detected by Spiesman (1995). However, 
they are very difficult to detect in the infrared because of the faint structure and the profusion 
of point sources. The data set used here is the high frequency component of the filtered signal 
shown in Figure 3. 

Here we use a method similar to that outline in the section for the smooth zodiacal background 
in order to co-add as many scans as possible to bring out the primary dustbands. The first three 
steps are the same as that for the smooth background. After obtaining the filtered scan, we use 
the scan for band 5 (12 microns) to fit polynomials to the ten-degree dustband pair. There are 
two different effect that have to be accounted for before co-addition. 

a) The distance between the pair (one above and the other below the ecliptic), denoted as (3, is 
plotted as a function of solar elongation. Increasing solar elongation results in lines-of-sight 
away from the sun towards the Asteroid Belt while smaller solar elongations look towards 
the sun. This introduces a parallax effect in the dust band pair such that the distance 
between the north and south bands increase with solar elongation. Therefore, the slope of 
the curve that defines the p with elongation is a measure of the parallax effect. This slope is 
now used to correct individual scans by using P at solar elongation of 90° as the 
normalizing value and is applied as ‘stretch’ or ‘shrink’ parameter for elongations less than 
and greater than 90°, respectively. 

b) The second correction is brought about by the fact that the dust bands are not parallel to the 
ecliptic. The inclination of plane of symmetry of the dust bands is brought about by long 
term planetary perturbations on the dust particles. The annual variation in the location of 
the dust bands is used to determine the inclination. Once the inclination is found, the 
correction is a ‘shift’ given by A Sin(X+<|)) where A is the amplitude, X is the ecliptic 
longitude of the Earth and <j> is the phase. The values at 12 microns have been obtained 
previously (Dermott 1992) and are as follows: 

Trailing: A = 1.25° ; <|> = 203.81° 

Leading: A = 1.14° ; <j> = 347.42° 
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Figures 8,9 andlO: The dust bands obtained from coadding the high frequency components at all 
elongations at three of the DIRBE near-IR bands. 


The results of the coaddition are shown in Figures 8, 9 and 10. It is very easy to detect the near- 
ecliptic dustband at all three wavebands while the ten-degree bands are visible only above or 
below the ecliptic at 1.2 and 2.2 microns. The zodiacal flux is a minimum near 3.4 microns 
hence it is far more difficult to detect them in band 3. The flux in the central dustbands are 
0. 1 1 3, 0. 1 1 5 and 0.054 MJy/Sr in bands 1 , 2 and 3 respectively. 

The data analyses presented in this report is part of an ongoing project to map the zodiacal 
cloud at all elongations and wavebands. The results will be published in Icarus in the near 
future. 

Work is currently being pursued to study the data set released by NASA after subtracting a 
zodiacal cloud model. The effectiveness of the model will be measured by comparing the 
results to that of the filtering process. In addition, results from the Midrange Space Explorer 
(MSX) at much smaller elongations (25-30°) are also being studied to sample the dust particles 
distribution very close to the sun. 
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